Visualizing paired electron pockets in the underdoped cuprates 
using ARPES in the presence of current 
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We propose an experiment to test a theory of the underdoped cuprates, which assumes that the pseudogap con- 
tains a phase-incoherent paired electron pocket and unpaired hole pockets. The proposed experiment involves 
angular-resolved photoemission spectroscopy (ARPES) measurements performed in the pseudogap regime and 
in the presence of a transport current running through the sample. The current leads to two main effects on the 
spectral function: First, even a weak current tilts the Fermi surface and is predicted to open up a part of the 
electron pocket if the energy of the incident photons is smaller but close to the pseudogap. Second, stronger 
currents suppress pairing of the electron pocket, which too can be observable by ARPES. The observation of 
these predicted phenomena, including their temperature and current dependencies, should clarify the central 
question about the existence of pairs in the enigmatic pseudogap region. 
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The effort to understand the behavior of copper oxide super- 
conductors [1] has lasted for decades. Theoretically, there are 
a large number of competing explanations [2] for the cuprates. 
They have been studied from the strongly correlated side at 
low doping or from weakly interacting models from the over- 
doped side of the phase diagram [3-14]. Despite serious ef- 
forts by theorists to explain the nature of the pseudogap, a 
consensus has not been reached and many key issues remain 
open. Among them is the question about the nature of the ex- 
citations in the pseudogap and in particular whether there exist 
preformed electronic pairs above the superconducting transi- 
tion. 

In this paper we propose and model an experiment that 
would help clarify this question and in particular may help 
to unambiguously determine whether the enigmatic pseudo- 
gap region contains "ghost" electron pockets. Specifically, 
we propose to perform angular-resolved photoemission spec- 
troscopy (ARPES) measurements in the presence of a current, 
which should produce useful data visualizing excitations that 
may otherwise be hidden due to pairing or gaps associated 
with other competing orders. 

While the usefulness of this proposed experiment is not lim- 
ited to a particular theoretical model, we focus on the specific 
model of paired electron pockets [15, 16], which we believe 
is strongly supported by three recent experimental develop- 
ments. First, recent work has observed a proximity-induced 
pseudogap [17] and supports the idea that the pseudogap is 
connected to paired quasiparticles. Second, the experimental 
discovery [18-20] of small Fermi pockets in the pseudogap 
phase of underdoped cuprates motivates a description which 
incorporates Fermi surface reconstruction. The apparent dis- 
crepancy between full Fermi pockets and the "Fermi arcs" ob- 
served by ARPES can be resolved with a model [21] consis- 
tent with what we consider here [16]. Third, there is a nodal- 
anti-nodal dichotomy observed in STM [22] and Raman [23] 
experiments. This is where nodal excitations have an energy 
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FIG. 1. (color online) A qualitative illustration of the experimen- 
tal setup and the predicted phenomenon: (a) the incident light and 
ejected electrons that initially leave the electron pockets hidden due 
to a pairing gap; (b) the electron-pockets are partially revealed in the 
ARPES signal due to a current running through the sample along the 
anti-nodal direction. The incident photon energy and direction of the 
current can be controlled to reveal hidden pairing. 



which decreases with decreasing doping, and anti-nodal exci- 
tations have a larger energy which increases with decreasing 
doping. 

These key ingredients, pairing, Fermi surface reconstruc- 
tion, and the nodal-anti-nodal dichotomy, are incorporated in 
a recent theoretical proposal describing the pseudogap phase 
of underdoped cuprates [16]. In this theory, the pseudogap is 
attributed to strongly-paired (but uncondensed) electron pock- 
ets in the anti-nodal region, mediated by orientational spin 
fluctuations (that suppress long-range Neel order in an ini- 
tially antiferromagnetic metal) and unpaired hole-pockets in 
the nodal regions (with a weak tendency for /?-wave pairing; 
the hole pairing is assumed to be destroyed by temperature 
in the pseudogap). This is shown to be equivalent to J- wave 
pairing when the Brillouin zone is unfolded [15, 16, 24]. Note 
that in the paired electron pocket scenario the key players are 
"ghost" electrons, which are always hidden from direct visual- 
ization due to a large pairing gap. Clearly, the main objective 
of an experiment that intends to test this picture is to unbind 
the electrons without destroying the underlying Fermi surface. 
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This situation might have been achieved in the quantum oscil- 
lation experiments, where the large applied field eff'ectively 
acts as a BCS upper critical field opening up the Fermi sur- 
face and yielding oscillations. However, ARPES measure- 
ments cannot be done in such a large field. Other possibilities 
include suppressing pairing by temperature and/or disorder, 
however they alone cannot reveal the "ghost" Fermi surface 
and will just smear it out. 

To test if the pseudogap is made of paired but uncondensed 
electronic excitations we propose an ARPES experiment in 
the presence of a current as shown in Fig. 1. There are two 
central reasons for applying a current. As we show, even 
a weak current will reveal the paired bands by shifting the 
spectrum. For larger currents complete depairing can occur 
opening up the "ghost" Fermi surface that should be visible 
by ARPES. This is pragmatic when the gap has already been 
reduced by raising the temperature so that the depairing cur- 
rent has a magnitude practical for experiment. 

Setup and model: We consider a hole-doped cuprate su- 
perconductor in the pseudogap regime subjected to a uniform 
current / as shown in Fig. 1(b). Our main goal is to determine 
key qualitative features of an ARPES spectrum measured in 
the presence of the current. We adopt a simplified model, 
which is consistent with the theory developed in Ref. 16 and 
contains the essential ingredients to reliably calculate the key 
features of the ARPES spectrum predicted by the theory. Ac- 
cording to this theory the pseudogap phase emerges from a 
strongly fluctuating critical antiferromagnetic state with re- 
constructed Fermi- surfaces consisting of electron-like pock- 
ets in the anti-nodal regions (nja, 0) and (0, nja) and hole-like 
pockets in the nodal regions (±7r/2a, ±7i/2a). The key formal 
ingredient is to describe the components of the physical elec- 
tron in a rotated reference frame set by the local spin-density 
wave (SDW) order (p = z^^o-o-o-'Zo-', where the bosonic spinon 
field Zo- defines the SU(2) rotation. This parametrization of 
the physical electron gives rise to an emergent gauge field, 
which plays a crucial role in the pairing of the fermions. The 
pseudogap phase is characterized by strongly ^--wave paired 
(but uncondensed) electron pockets and hole pockets [25] that 
remain unpaired (a weak /?-wave pairing in the hole pock- 
ets is assumed to be suppressed by temperature). Here, we 
adopt a simple minimal model for the pseudogap phase which 
nevertheless captures the key features necessary to address an 
ARPES experiment in the presence of current: 

-Y,Wlfi,_+h.c.] + ^, (1) 

k 

where fk± and hk± are the annihilation operators for the 
electron-like and hole-like excitations with momentum k, re- 
spectively, and ± labels the charge associated with the emer- 
gent gauge field. The sum over k is performed over the re- 
duced Brillouin zone. We note that fka becomes the physical 
electron in the SDW ordered state, where the spinon field Zo- 



condenses and the ± indices become equivalent to the elec- 
tron spin indices. Pairing in the electron pockets is included 
at the mean-field level by introducing a real ^-wave pair po- 
tential A. The spectra are given by ^kf,h = i^k + ^k+Q)/2 ± 
[(^k - ^k+af + ^<P^] ^^^/2, where = -2^(cos k^a + cos kya) - 
4f cos kxa cos kya - ji. Here, Q = (n/a, n/a), and we will take 
f = -03t and ji = -0.6t [26]. The quantity (f is the uni- 
form SDW order parameter, but we stress that cp is used here 
merely to parameterize the underlying Brillouin zone folding, 
and that in reality, spin fluctuations are expected to suppress 
long-range antiferromagnetic order in the pseudogap phase. 
In the above, A> Ois the efl'ective attractive interaction for the 
electrons generated by the gauge fluctuations. Since the pseu- 
dogap phase is not characterized by long-range SDW order 
the attractive interaction mediated by the gauge fluctuations 
is, in principle, long-ranged [16]. 

We emphasize here that a calculation of the current within 
model (1) would give rise to a phase-coherent superflow, 
which is not correct. In the pseudogap, pairs have formed but 
transport is in terms of holes and phase-incoherent pairs of 
charge 2e. However, the excitation spectrum in the presence 
of current should be correctly obtained from these expressions 
and can reliably calculate the spectral function, which is the 
central quantity of interest. 

Spectral function for a current-carrying pseudogap: We 
now study modifications to the ARPES spectrum when a cur- 
rent is applied to a cuprate superconductor in the pseudogap 
regime. Although ARPES is a probe of the physical electron 
spectral function we will compute the spectral functions for 
the electron- and hole-like excitations in this work. In the 
presence of spinon fluctuations we cannot, in principle, expect 
the physical electronic spectrum to be identical to that of the 
electron-like excitations [21]. However, as we will show, the 
application of even a weak current can cause sections of the 
electron pockets to appear or disappear. Furthermore, a strong 
enough current will completely depair the electron pockets 
and reveal their entire Fermi surfaces. Despite our simplified 
treatment, these key modifications should persist also in the 
ARPES spectrum as well. 

The current endows the Cooper pairs with a center-of-mass 
momentum qj, and the single-particle Green function in the 
presence of a uniform current can be written as [27, 28] 

Gkiicon) = ^{icon " Vk ' ^y)! - hfh + Afi] ^ , (2) 

where cOn = (2n -h \)7ikBT are the fermionic Matsubara fre- 
quencies, Vk = d^kf/hdk, checks indicate matrices in Nambu 
space, and f / are the Pauli matrices. The gap is then obtained 
from the following self-consistent condition. 

Ay l-nF{El)-nF{El) 

N ^ 2Ek ' ^ ^ 

where nf is the Fermi-Dirac distribution function and = 
Ek ± fiVk • qj = + A^]^/^ ± fiVk ' qj- The retarded Green 
function for the -i- electrons, G^^{(jS), then defines the corre- 
sponding spectral function A (a;) = -ImG^^(^^)/7r. One then 
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finds, 

Ak+(oj) = s\6(n(jj - ^kh) 

+ rl [ulSihoj - El) + vldifioj + ED] , (4) 

where = [1 + (^^ and ^ = [1 - - 

^k+Q)/('^kf - and the Bogohubov coherence factors 

are given by = (1 + f^^/Ek)/2 and = (1 - f^j^/Ek)/2. 
Ahhough, interaction and disorder eff'ects beyond the scope 
of our simpHfied model will inevitably give rise to non-trivial 
self-energy corrections and modifications to (4), we are even- 
tually interested in looking at the robust equal-energy contours 
of the dispersions rather than the level widths, which are sen- 
sitive to a self-energy. Therefore, (4) should correctly convey 
the main qualitative eff'ects of the applied current on the spec- 
tral function. 

Depairing due to current. The gap solution to (3) as a func- 
tion of qj along the anti-nodal direction is shown in Fig. 2. 
As illustrated in the inset, the gap is plotted for the spectrum 
^kf in (a) and, for comparison, for the regular parabolic spec- 
trum ^ = ff-k^/2m - fi in two dimensions in (b). The gap is 
normalized by Aq, which is its value at zero temperature and 
zero current. The values (f = 03t and t/A = 0.7 were used in 
(a), and A^o^ ~ 0.34, where A^o is the density of states at the 
Fermi level, was used in (b). 

We find that the Fermi velocity for the electron pockets 
is hvF/a ^ 1.3^. The depairing scale should then be set by 
q^j ~ Ao/hvF- As we see in Fig. 2(a), at zero temperature, the 
gap remains robust up to this scale and then shows a steep but 
gradual decrease. A general rounding of the features is ob- 
served at finite temperatures. For the parabolic dispersion at 
r = 0, the gap remains robust for hvFqj < Aq, then sharply 
drops to zero at the critical depairing scale hvFq^j ~ Aq. At 
finite temperatures, the gap monotonically and smoothly goes 
to zero. We note that the shoulder for the T = result 
in Fig. 2(a) appears because the gap equation (3) sums the 
doppler shift over a non-trivial Fermi surface corresponding to 
the dispersion ^kf- The depairing scale Aq is set by the pseu- 
dogap scale, and can be reduced by raising the temperature 
of the sample until Ao(r) is sufficiently small. Here, Ao(r) 
is the gap at temperature T and zero current. To the best of 
our knowledge, depairing eff'ects in two-dimensional conven- 
tional superconductors have not been thoroughly investigated, 
since the initial interest in the early days of BCS theory was in 
three-dimensional materials. Many investigations of three di- 
mensional ^--wave [29-31], two dimensional J-wave [32, 33] 
and other exotic pairing [34, 35], however, have been done. 

For a weak current, one can obtain an approximate analytic 
expression for the gap. Solving (3) perturbatively, one finds 
that the gap as a function of current (in the anti-nodal direc- 
tion) is given by: 

A(qj) ^ Ao(T) [l - rifiVFqj/ Aof] , (5) 
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FIG. 2. Solution to the gap equation (3) as a function of along the 
anti-nodal direction. In (a), the solution is obtained for electrons with 
spectrum ^kf (with the Fermi surface shown in the inset), while in (b) 
we consider the regular isotropic spectrum = ff-k^ llm-ji (with the 
circular Fermi surface as shown in the inset). The gap is normalized 
with respect to its value at zero temperature and zero current. 



with Xk = ElllT and El = [f^^ -h Al(T)Y^^. For the T = 
0.45 Ao resuk in Fig. 2(a), we find Ao(r = 0.45 Aq) ~ 0.7 IAq. 
In this case, y ^ 0.58, and the perturbative treatment here 
describes the initial downturn in the gap. 

Spectral function results: From (4) we see that the spectral 
function is non-zero when either (i) hoj - hvk • qj - Ek = 
or (ii) fioj - fiVk ' qj + Ek = 0. We see that the eff'ect of the 
Doppler shift is to move the energy scale fioj away or towards 
the quasiparticle bands ±Ek depending on the sign of fiVk • qj. 
While the depairing current scale set by q^j ~ Ao/hvF can 
require a large current, all that is required here is to use a small 
^Vk • qj to just cross into either of the ±Ek bands. In Fig. 3, 
we plot the spectral function (4) at photon energies just below 
[Fig. 3(a)] and above [Fig. 3(b)] the upper paired band for zero 
and finite current. We consider the case where the current 
is applied in the anti-nodal direction, and also in the nodal 
direction [Fig. 3(c)]. Plotted on the left are a cut of the paired 
band dispersion along the = line in the Brillouin zone 
(indicated by the red dashed line on the right in (a)), and the 
photon energies are indicated by the dotted lines. In (a), the 
photon energy is set inside the gap at zero current and hence 
the electron pockets are initially not observed. However, once 
the current is applied, the paired spectrum at k points where 
Vkis parallel to qj is shifted up in energy while the spectrum 
at k points where Vk is anti-parallel to qj is shifted down in 
energy. This leads to the appearance of some sections of the 
electron pockets. In Figs. 3(b),(c) the photon energy crosses 
both the hole and the upper paired bands. The application of 
the current here leads to the disappearance of some sections 
of the electron pockets. 

We remind the reader that the spectral functions plotted in 
Fig. 3 (see (4)) are for the electron-/;!^ excitations residing in 
the anti-nodal regions of the Brillouin zone (we are ignoring 
the hole-like spectral function here). For the ARPES spec- 
trum, the spectral function must be computed in the physi- 
cal electron basis obtained from rotating back our electron- 
like excitations using the spinon field Zo-- In the absence of 
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FIG. 3. (color online) Spectral function at zero and finite current for 
photon energy hoj just below [see (a)] and above Ek [see (b)]. As 
indicated by the white arrows, the current is applied along the anti- 
nodal direction in (a) and (b), and the nodal direction in (c). Paired 
bands on the left are plotted along the ky = line in the Brillouin 
zone, as indicated by the red dashed line. In principle, an arbitrarily 
small q j can reveal or erase parts of the electron pockets as long as 
the photon energy is tuned close to the edge of a paired band (see text 
for details). 

long-range Neel order, these spinons are uncondensed and the 
spinon fluctuations are known to play an important role in the 
ARPES spectrum [21]. Nevertheless, the appearance and dis- 
appearance of the electron-like Fermi surfaces found in this 
work should have a clear signature in the ARPES spectrum 
as well. This is because to lowest order the physical electron 
spectral function is simply a convolution of the spectral func- 
tions of its composite electron-like and spinon particles [21]. 
The appearance or disappearance of the Fermi surface for the 
fermion component should then directly afl'ect the convolved 
spectral function. 

Heuristically, we now provide the relationship between the 
drift momentum qj and the current density based on the 
Drude model. The total current J = Jb + Jh '^^ composed 
of holes from the hole pockets Jh = = (nhe^Th/mh)E 
and from uncondensed pairs Jh = (TbE = (niy(2e)^Tiy/2me)E. 
Here, rih is the number density of uncondensed pairs and rih 
is the density in the hole pockets. In the pseudogap phase we 
have a mix of two dissipative fluids whose conductivities will 
add J = {(Tb + ChjE. At steady state, qj = leEzb. This 
means that for a given / the drift momentum is related to total 
current via qj = 2ejTh/(crh + ct/^). 

The simplest estimate for the critical depairing current den- 



sity is Jc = 2enAo/mVf, where n is the density of charge 
carriers. Specific tight-binding models for Bi2Sr2CaCu208 
(BSCCO) can be taken from Refs. 26 and 36. In BSCCO, 
vf ~ 250km/s [37], and the eflTective mass is ^ 3-4me [38]. 
A typical carrier density in BSCCO [39] is ^ lO^^cm"^. 
Therefore, Jc ~ 65AoA/cm^ for Aq in units of Kelvin. The 
energy scale for Aq can be estimated as follows. The ARPES 
pseudogap can be identified by two energy scales [40], with 
a low-energy scale given by the location of the leading-edge 
midpoint and a higher-energy one by the position of a broad 
peak. The broad maximum of the ARPES spectra are between 
100-200meV [41] and the leading-edge midpoint are between 
around 20-30meV [40]. Scanning tunneling spectroscopy [42] 
in BSCCO also sees a pseudogap feature around 70-90meV. 
Ref. 22 also sees a high energy feature in their data in the 
60-lOOmeV range. In Ref. 16, the pairing scale is identi- 
fied as Ao ~ 200K consistent with the low energy features 
in ARPES and the downturn in uniform spin susceptibility 
occurs around 300-400K. The observation of a proximity in- 
duced pseudogap [17] sees a feature in the normal density of 
states at around sample bias lOmeV. This would correspond 
to a gap smaller than the low energy feature, but due to gap 
attenuation in the proximity eflTect, this might be consistent 
with the small pseudogap scale. Using Aq ~ 200K, one finds 
/c ~ 1.3 X 10"^ A/cm^. The magnitude of this estimate is not 
unreasonable if we note that a bulk critical current density of 
2-5x10^ A/cm^ has been obtained for BSCCO [43, 44]. 

In conclusion, we have shown that the application of a weak 
current during ARPES in the pseudogap state provides a test 
of the model proposed in Ref. 16, by either revealing sections 
of the electron pockets, or causing them to disappear, depend- 
ing on the incident photon energy and current direction. We 
have also solved the two dimensional gap equation in the pres- 
ence of a current, establishing the critical depairing current 
scale hvFq^j = Aq for a circular Fermi surface but depends 
on Fermi surface geometry in general (Fig. 2). To pragmati- 
cally reach the depairing current, Aq needs to be reduced by, 
for instance, temperature. While the photo-ejected electron's 
trajectory will be modified by the presence of a weak field in- 
duced by the current, it can be accounted for. We have thereby 
provided a falsifiable test for a model of the cuprate pseudo- 
gap. 
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